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^ Abstract 

P5 , Contact P-wave interactions connected to the Larmor interaction of a magnetic 
dipole and Thomas spin precession in the filed of an electric quadrupole are de- 

ff^ ■ scribed and their implications for spectroscopy of exotic ri~-atoms are studied. In 

i order to evaluate the magnitude of the contact P-wave interactions as compared to 

I the conventional long-range interactions and the sensitivity of spectroscopic data 

■ to the J7~-hyperon quadrupole moment, we consider 2P states of $7" atoms formed 

m 



with light stable nuclei with spins / > 1/2 and atomic numbers Z < 10. The energy 
I level splitting caused by the contact interactions is 2 — 5 orders of magnitude smaller 

■ than the conventional long-range interactions. Strong decay widths of atoms 

(3 ■ due to reactions pQ~ AH" and pQ~ SH, induced by t-channel kaon exchanges, 

are calculated, fi" atoms formed with the light nuclei have strong widths 5 — 6 
orders of magnitude higher than splitting caused by the contact interactions. The 
low-L pattern in the energy spectra of intermediate- and high-Z atoms thus 
5^ \ cannot be observed. The quadrupole moment can be measured by observing X- 

rays from circular transitions between high-L levels in fi" exotic atoms. The effect 
of strong interactions in ^'^^Pbfi" atoms is negligible starting from L ~ 10. The 
contact P-wave interactions exist in ordinary atoms and /x-meson atoms. 
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1 Introduction 



In the early 70's, Goldhaber and Sternheimer [T] proposed to measure the 
i7^-hyperon magnetic and quadrupole moments by detecting X -rays from 
circular transitions of ^"-hyperons captured to atomic orbits. Such method 
has been successful in extracting experimentally the magnetic moment of the 
E~-hyperon from fine splitting in S~ exotic atoms |2j. The fl~ magnetic mo- 
ment has been measured with high precision, however, by other techniques 
[311]. The measurement of the fl~ electric quadrupole and magnetic octupole 
moments remains an open problem. 

Among the decuplet baryons, the ^7^-hyperon has weak decays only and a 
small width. It appears to be a suitable candidate for measurement of the 
static quadrupole moment. The measurement of the f2~-hyperon quadrupole 
moment would be helpful to understand better hadron structure and proper- 
ties of quark interactions. 

The fl~ exotic atoms are discussed in Refs. OHEISE]- The f2~-hyperons are 
produced experimentally as relativistic particles. Stopping Q~ is a hard exper- 
imental task, since ^"-hyperons in matter dissolve to lighter hyperons by ex- 
changing ii"-mesons with surrounding nuclei. The reaction K~p K~K^Q~ 
at threshold is in particular not well suited for producing slow ^7^-hyperons 
[To] . During the last three decades, there has been no progress in experimental 
studies of Q~ exotic atoms. 

Two events of stopped S~-hyperons in light emulsion nuclei at KEK have been 
interpreted as S~ atomic states bound with ^^C [TT|IT2] . Future experiments 
for producing high rates H- hyperons at GSI are discussed in Refs. [TSfTl] . 
Properties of E~ atoms are discussed in Refs. [TUfT5] . 

Recently Karl and Novikov [T6p7] made an interesting observation on the 
existence of a contact P-wave interaction of two quadrupoles and proposed to 
measure the f2~-hyperon quadrupole moment from the hyperfine splitting of 
P-wave Q~ atomic states. The f2~-hyperon is the only (almost) stable particle 
which can form bound states with a high-spin nucleus to exhibit quadrupole- 
quadrupole interactions. 

The fine and hyperfine splittings in atoms relative to the ground state energy 
are of the order (f/c)^ ~ {ctZY, the hyperfine interaction is suppressed addi- 
tionally by a factor ~ rrie/M, where rrie is the electron mass and M is the mass 
of nucleus, the Lamb shift is of the order ~ a{aZY log while the contact 
P-wave interaction is of the order a{aZ)^. The additional smallness ~ a'^Z as 
compared to the dominant terms might be compensated in individual cases by 
large quadrupole moment of a high-Z nucleus and/or a specific pattern of the 
quadrupole-quadrupole splitting. In this work, we analyze hyperfine splitting 
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in Q,~ atoms by comparing numerically the magnitude of various interactions 
in Q~ atoms, formed with light stable nuclei, including spin-orbit interactions, 
spin-spin interaction, quadrupole-orbit interactions, which are of order (f/c)^, 
and contact P-wave interactions of order (u/c)'*. 

The measurement of energy splitting is possible provided widths of the corre- 
sponding energy levels are small. We calculate the strong decay widths of pQ~ 
exotic atoms with arbitrary principal and orbital quantum numbers and give 
rough estimates of the strong decay widths of fl~ exotic atoms formed with 
high-Z nuclei. 

The outline of the paper is as follows: In the next Sect., we discuss config- 
uration mixing and exchange current contributions to quadrupole moments 
of the decuplet baryons and other static observables of baryons. In Sect. 3, 
a description of various interaction terms in bound systems, which appear in 
the nonrelativistic expansion of the one-photon exchange interaction potential 
between two high-spin particles, is given. The isotope dependence of the spin- 
orbit interaction is discussed. A contact P-wave electric quadrupole - magnetic 
dipole interaction is described and its magnitude is estimated and compared 
to other interactions in atoms. In Sect. 4, we describe the calculation of 
strong decay widths of pQ~ exotic atoms due to the processes pQ~ —>■ AS, SS. 
Strong decay widths of 2P states of fl~ exotic atoms with light nuclei up to 
^^F are found to be up to three orders of magnitude higher than the dominant 
long-range interactions. Estimates made for circular transitions in ^°^Pbf2~ 
exotic atoms give small strong decay widths starting from L = ri — 1 ~ 10. 

In Conclusion, we summarize the results. 



2 Configuration mixing vs two-body exchange currents 



Quark models are known to be very successful in the description of hadron 
properties. The one-gluon exchange describes the quadrupole moments of the 
decuplet baryons |18fl9f20ll21f22f23ll24j and the non-vanishing neutron charge 
radius |24f25f26ll27f28f22j . In the framework of the Isgur-Karl nonrelativistic 
quark model, these quantities are simply related [23] : 

2 

Qa+ = 3r„ \cM , (2.1) 
5 

where Qa+ is the A'^'-isobar quadrupole moment and r„ is the neutron charge 
radius, determined by configuration mixing (CM) in the baryon wave functions 

^ The quadrupole moments of nuclei increase with Z roughly as Z^/'^, so the contact 
P-wave quadrupole-quadrupole interaction is well approximated as (v/c)^ ~ (aZ)^. 
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as illustrated on Fig. 1 (a). 



Let us discuss the status of CM effects in terms of the v/c expansion. Spin-spin 
forces in the Fermi-Breit potential are of the order Vss ~ l/(m^r'^), where m 
is the constituent quark mass. The corresponding perturbation of the baryon 
wave functions is of the order 5'^ ~ ~ / {m?uor'^), where A£' ~ (cu is 

the oscillator frequency). Thus the neutron charge radius and the quadrupole 
moments are of the order l/[rn?ujr). The ratio between the neutron charge 
radius ~ 1/ {m^ujr) and the proton charge radius ~ 1/ {muj) becomes 
''^nl^p ~ ^/mu/m ~ t>/c, where we have used the relations ~ muj for 
an oscillator and •pjm ^ v j c. CM effects in the quadrupole moments and the 
neutron charge radius are therefore of the order vjc. Refs. p!8lll9|20|23|22j and 
Refs. [25|26f22j provide the calculations of Qa+ and r^, respectively, using the 
nonrelativistic quark model and Ref. [27] provides the calculation of r\ using 
MIT bag model. Refs. [TMM2nll23ll25ll26ll?^ evaluate the CM. 

Two-body exchange currents (EC) in bound systems contribute to observ- 
ables also. They are associated to tree level Z-diagrams of the noncovariant 
perturbation theory, shown on Fig. 1 (b). 




(a) (b) 

Fig. 1. Configuration mixing diagrams (a) and exchange current diagrams (b) con- 
tributing to an observable marked by the crosses. The solid lines are quarks, the 
dashed Unes are gluons and mesons. 

EC corrections to the charge density operator are of the order 1 / (m^r^) [29)130] , 
so the corresponding corrections to the quadrupole moments and the neutron 
charge radius ~ l/(m^r). The ratio between EC correction to the neutron 
charge radius and the proton charge radius becomes ~ ujm/ {m?r) ~ (f/c)^. 
EC corrections to the charge density operator and therefore to the quadrupole 
moments and the neutron charge radius are of the order (f/c)^. In the frame- 
work of the nonrelativistic quark model, one can expect that EC effects are 
small as compared to CM effects for observables related to the charge density 
operator. 

Precise measurements of the transition quadrupole moment A+ p'-f give 

^ A subclass of Z-diagrams shown on Fig. 2(d) of Ref. [27j vanishes. The remaining 
class of diagrams shown on Fig 2(b,c) corresponds to configuration mixing. 
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a value Qa+p7 = -0.108 ± 0.009 ± 0.034 fm^ [3T] significantly higher than 
values predicted by the nonrelativistic quark models p^HTO] based on evalua- 
tion of CM alone with realistic quark core radii. One can expect that static 
quadrupole moments are undervalued too. Buchmann, Hernandez and Faessler 
[32] conjectured that EC effects in observables related to the charge density 
operator are dominant. If one neglects CM effects and keep EC effects, one 
gets relation [32] 

Qa+ = rl \ec , (2.2) 
which gives a higher value for the quadrupole moment of the A. 

Relativistic quark models sum up the v/c series. It is thus instructive to com- 
pare the nonrelativistic quark model predictions with relativistic models. The 
experimental value of Q^+p-y appears to be three times higher than prediction 
of Ref. [2T] based on the chiral bag model with account taken of CM and EC 
effects Li|. The MIT bag model calculation of Close and Horgan [27] where CM 
effects are included only gives the neutron charge radius much smaller than 
that obtained in Ref. [2S|- This result agrees qualitatively with the conjecture 
of Buchmann, Hernandez and Faessler [52] on the dominance of higher order 
v/c terms in observables related to the charge density operator. In the chiral 
bag model, Q^+p-y and are still undervalued. EC corrections to Q^+ and 
are calculated Refs. [32p^ using the nonrelativistic quark model and in Refs. 
[2Tp8f33j using MIT and chiral bag models. 

It is known that one-gluon exchange contributes to magnetic moments of 
baryons. EC contributions to the current density operator of nonrelativistic 
systems can be obtained from the Fermi-Breit potential by the minimal sub- 
stitution p — > p — eA and taking derivative of the potential over A. Magnetic 
moments of composite systems receive corrections ~ l/(mr) ~ f/c. The 
corresponding CM corrections due to the orthogonality of the space part of 
the quark wave functions are proportional to Sjj/n ~ (^^)^ ~ {v/cY. In 
the framework of the nonrelativistic quark model, one can expect that EC 
corrections are large as compared to CM corrections when observables are re- 
lated to the current density operator. Such a premise does not contradict to 
observations. CM corrections to baryon magnetic moments are calculated in 
Refs. [Mf2^ in the nonrelativistic potential model and in Refs. [35|37f38f33f39] 
using the MIT bag model. EC corrections to baryon magnetic moments are 
calculated in Refs. [36|32II24] using the nonrelativistic potential model and in 
Refs. [371l38ll33|l39] using MIT and chiral bag models. 



^ Ref. [21] treated gluon and quark fields classically. As shown in Ref. [33j, sum- 
mation of tree diagrams of the perturbation theory is equivalent to solving the 
classical equations of motion for gluon and quark fields, Z-diagrams are contained 
in the lower components of Dirac bispinors describing interacting quarks. Owing to 
quark self- interactions, results of Ref. [21] give quantum predictions. 
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Modern relativistic quark models [lO] treat exchange effects between quarks 
on the basis of the covariant perturbation theory without exphcit separation 
to CM and EC effects. 

The mesons exchange effects are suppressed by the mass of the exchanged 
mesons. One can expect that in the f2~-hyperon the role of mesons is less 
important. The measurement of the f2~-hyperon quadrupole moment can be 
helpful to differentiate the gluon and meson exchange effects. 



3 Fine and hyperfine interactions 



Fine and hyperfine interactions in ordinary atoms are described in standard 
textbooks (see e.g. [^T]H2] and others), while specific features of fl~ exotic 
atoms are discussed in Refs. [5H6H71l8ll9] . The discussion is, however, restricted 
to spin-zero nuclei. For a high-spin nucleus, the pattern of energy levels is 
more complex due to the presence of higher nuclear multipole moments. 

In this section, we summarize the known facts about interactions of high- 
spin particles, contributing to the energy level splitting to order (f/c)^, and 
describe contact P-wave interactions of order (f/c)^. Our purpose is to check 
the numerical magnitudes of various contributions to the energy splitting of 
an Q~ atom in L = 1 states with an intermediate mass nucleus. In Sect. 4, 
we compare the energy level splitting with widths of Q~ exotic atom due to 
reactions pQ~ — > AS, ES. 



3.1 Isotope effect in spin-orbit interactions 



The binding energy of electrons in atoms depends on the mass of nuclei, M, 
through the reduced electron mass 



In the hydrogen-like atoms. 



mM , , 

^ = T7- 3.1 



En = (3.2) 



where n is the principal quantum number. The isotope effect in the energy 
levels of hydrogen-like atoms is measured experimentally and described in 
standard textbooks (see e.g. [12]). 
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The spin-orbit splitting in atoms represents a relativistic effect ~ (f/c)^. In 
the hydrogen-hke atoms, 

^ = ^. (3.3) 

c n 

Corrections ~ me/M ~ 5 ■ 10~^ to the spin-orbit sphtting in ordinary atoms 
are usually not discussed. However, exotic atoms such as antiproton and S~- 
hyperon atoms are created and studied in the laboratory. In Ref. [2], the S~ 
magnetic moment has been measured from fine structure in S~ exotic atoms. 
In exotic atoms, the isotope effect becomes important. 

The LS potential consists of two parts. The Larmor part is connected to the 
interaction of the magnetic moment of the bound particle with the magnetic 
field generated by the nucleus in the co-moving frame of the bound particle. 
The second contribution is related to the Thomas spin precession. We thus 
write 

Uls = Uts + Uls. (3.4) 

The electrostatic potential created by a nucleus at rest with charge —eZ has 
the form 

eZ 

$ = , 3.5 

r 

where e = — |e| is the electron charge and r = |x|. The electric field equals 

E = -V$. (3.6) 

In the co-moving frame of the bound particle, the magnetic field can be found 
using the Lorentz transformation: 

B' = -V X E = V^L, (3-7) 

m r dr 

where v = vi — V2 = pi/m — P2/M = p/m' is the relative velocity, p = pi = 
— P2 in the center-of-mass frame, L = x x p where x = xi — X2. The indices 
1 and 2 refer to bound particle and nucleus, respectively. 

The Larmor component of the spin-orbit interaction potential becomes 

£s ■ B' = 

S 2mm'r^ 



UEs = . B' = -^L . S, (3.8) 



where 

is the magnetic moment and g is the gyromagnetic ratio of the particle (for 
electron S = 1/2 and g = 2). 

The angular frequency of Thomas precession, fix, is related to the angular 



7 



frequency of the orbital motion, uj: 



Or = a;(l-7), (3.10) 

where 7 is the Lorentz factor of the moving particle. Equation fl3.10p is derived 
in Appendix A. The Hamiltonian producing the spin precession fl3.10p is given 

by 

Ht = ^t- S. (3.11) 
To the first order in (f/c)^, one gets 

Uls = nr-S = -^l..S. (3.12) 

Here, one used L = xi x pi + X2 x p2 = (xi — X2) x p = mrRuj, where 
R = m'r/m is the distance from the center-of-mass of the system to the 
particle, p = mRu, and mv'^/R = aZjr^ for particle on a circular orbit. In 
an external potential of scalar type with respect to the Lorentz group, the 
probing particle experiences the Thomas precession only, so its spin-orbit po- 
tential takes the form of Eq. fl3.12p . The isotope dependence of scalar-exchange 
potentials, implied by Eq. fl3.12p . is in agreement with Ref. [13]. 

The spin-orbit potential in the Coulomb field takes the form 

f/.. = ^(^-i)L.S. (3.13) 

The Fermi-Breit potential for the f2~-hyperon atom of Ref. [7] contains Eq. fl3.13p 
A distinct isotope dependence of the spin-orbit potential is used in Refs. [T1IT7] . 
For ^ = 1/2 and ^ = 2 Eq. (l333|) is in agreement with Ref. 01], Chap. IX. 

For high-Z atoms, the Dirac equation is usually used, modified to include the 
anomalous magnetic moment of the particle and with m replaced by m' . Borie 
[15] developed efficient numerical schemes for calculation of energy eigenstates 
of relativistic atoms including the nucleus recoil corrections. The spin-orbit in- 
teraction obtained in Ref. [IS] by the non-relativistic reduction of the modified 
Dirac equation is in agreement with Eq. fl3.13p to order 1/A. This accuracy is 
sufficient for extracting the S~ magnetic moment from the spin-orbit split- 
ting of the high-Z exotic atoms Pb-S^ and W-E~ [2]. In low-Z atoms such as 
p-^He or K^-p the corrections are important. 

The magnetic field created by the orbital motion of the particle acts on the 
magnetic moment of the nucleus. The nuclear spin I experiences the Thomas 
precession as well. The spin-orbit interaction potential has the form: 



The magnetic moment of the nucleus is defined by 



^^^^ T 

For high-Z nuclei, the Larmor contribution to Ujl is of order of unity, whereas 
the Thomas precession is suppressed as 1/A. In the potential Uls these con- 
tributions are of the same order in 1/A. 



3.2 Spin-spin interaction 

The long-range interaction of two magnetic dipoles is well known (see e.g. [41], 
Chap. XVI): 

The irreducible tensor r"^(a, b) with space indices a, /3 = 1, 2, 3 is defined by 
r"^(a, b) = a"6^ + a'^fo" - ^a • b(5"^. (3.17) 
Its properties are described in Appendix B. 

3.3 Quadrupole- orbit interactions 



The electrostatic potential created by a spin-/ nucleus gives rise to hyper- 
fine splitting connected to the interaction of the nucleus electric quadrupole 
moment with electrostatic field created by the orbital motion of the bound 
particle. In terms of the nucleus spin I, the nucleus electric quadrupole mo- 
ment has the form 

Qz' = - J^^^\f '{lJ)- (3.18) 

It is normalized by 

< n\Q'i\n >= -eQz- (3.19) 
For nuclear electric quadrupole moments, experiments quote Qz i.e. the proton 
charge, — e, is usually omitted. 

The quadrupole-orbit interaction has the form (see e.g. [42] ) 

The hyperfine splitting caused by the potential Uq^l is used to measure the 
electric quadrupole moments of nuclei from spectroscopy of ordinary atoms 
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[12] and X-ray spectroscopy of vr-, K- and /i-meson exotic atoms (see e.g. 

The nucleus Coulomb field interacts with the quadrupole moment of VL~ 

«"=2S(I|3T)-"(S.S)- (3-21) 
It is normalized as follows: 

< SS\Q''\SS >= eQ. (3.22) 



The interaction potential has the form of Eg. (13.201) with the replacements 
a ^ aZ, Qz ^ Q, and I S: 



3.4 Quadrupole- spin contact P-wave interactions 



Contact S'-wave interactions of baryons originating from meson exchanges are 
usually omitted from the start, since those interactions are made to vanish by 
the repulsive core anyway. The S'-wave contact terms, generated by photon 
exchanges, can be set equal to zero in the baryon-baryon potentials also. 

Contact P-wave interactions depend on the gradient of wave function at the 
origin. Assuming the wave function is suppressed as ^'(O) ~ exp(— A/2m?7o6) -C 
1 where Uq is height of the core and b is the core radius, one finds that \l/(0)' is 
suppressed as ^'(O)' ~ ■\/2mUoaB exp{—^y2mUob), where is the Bohr radius, 
qb ^ h. One sees that suppression of the contact P-wave interactions is less 
effective and can in principle be compensated by the large factor y/2mUoaB- 

We analyse contact P-wave interactions quantitatively and, as a first approx- 
imation, assuming no effect from the repulsive core exists at all. 

Let us consider the quadrupole part of the electrostatic field 

6 r 

Tensor V'V^V^^ entering this expression can be split into two parts with 
angular momenta L = 3 and L = 1: 



V"V^V^- = T°f^ + T['^f^, (3.25) 
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where 

5 7* 5 



In the co-moving coordinate system of the f2 hyperon the induced magnetic 
field has the form 

B'^ = -e^'^^v^'E^ (3.26) 

where v'^ = p'^/m' is velocity, p'^ momentum and m' reduced mass of Q~, e^'^"' 
is the totally antisymmetric tensor such that e^^^ = 1. The Larmor interaction 
energy of the VL~ magnetic moment with the magnetic field is given by 

U^s = (3.27) 
The contact part of the interaction has the form 

^£5 = ^^e'^^QfS\p^V-6{^)) + (V"<5(x))p^). (3.28) 
The term entering the brackets, being averaged over the L = 1 state, gives 



/ d^Y*^,{n)R,,ir) (p"(V^5(x)) + (^^^(x))^") 

where Rj^i^{0) is the derivative of the radial wave function at the origin Eq.( ]B.13|) . 
We thus obtain 

U'J^s = T"/^(I,I)r"^(S,L)i?:,^,(0). (3.29) 

40mm' 2/(2/ - 1) V . ; V . ; mv ; v ; 



For evaluation of the Thomas precession component of the interaction energy, 
we determine from Eq.(A.22) the Thomas precession frequency, 

riT ^ -^vi X vi, (3.30) 

and use Eq. 03.111) . For circular motion is in agreement with the nonrela- 
tivistic limit of Eq. fl3.10p . Taking for E expression fl3.24p and substituting its 
contact part into the equation of motion md\i/dt = eE, and further, d'Vi/dt 
into Eqs. (13.301) . we obtain the Thomas component, Ucj^si contact in- 

teraction potential. If we would use for E expression (III. 6), we could get 
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Table 1 

The magnitudes of long-range and contact P-wave interactions and strong decay 
widths of atoms in 2P states, formed with light stable spin I > \ nuclei. The 
experimental values of the nuclear dipole magnetic moments and the nuclear 
electric quadrupole moments Qz are taken from Ref. [51], errors are not displayed. 
The multipole moments of O": ^ = -2.02 n.m. [HI], Q = -2.8 x IQ-^ fm^ [23j. C/i^L 
is the maximum value over F of the root mean square of eigenvalues of the matrix 
elements < FJ|C/[^1|FJ' >, Umax is defined similarly. For ^H, the strong decay width 
is given by Eq. (|4.59p . in other cases T are estimates based on Eq. (j4.63p . 



Nuclei 




2H 


3H 


^He 


6Li 


^Li 


9Be 


1 


1 /2 


1 






1 


3/9 


3/9 


jJLz [n.m.] 


2.79 


0.86 


2.98 


-2.13 


0.82 


3.26 


-1.18 


Qz [fm2] 




0.29 






-0.08 


-4.06 


5.3 


f/rk'ix [keV] 


9 X 10-6 


3 X lO-'' 


10-4 


6 X 10-2 


2 X 10-2 


0.2 


1 


t/,^L [keV] 


4 X 10-10 


7 X 10-9 


4 X 10-9 


10-7 


2 X 10-6 


3 X 10-5 


3 X 10-4 


r [keV] 


6.3 X 10-6 








10 


20 


100 


Nuclei 


lOg 


iiB 


13c 


14n 


15N 


1^0 


19p 


/ 


3 


3/2 


1/2 


1 


1/2 


5/2 


1/2 


HZ [n.m.] 


1.80 


2.69 


0.70 


0.40 


-0.28 


-1.89 


2.63 


Qz [fm2] 


8.47 


4.07 




2.00 




-2.58 




f/rk'ix [keV] 


2 


2 


0.5 


3 


1 


3 


2 


t/,^L [keV] 


4 X 10-2 


8 X 10-4 


3 X 10-^ 


6 X 10-2 


3 X 10-^ 


2 X 10-2 


10-2 


r [keV] 


350 


500 


10^ 


2 X 102 


4 X 102 


7 X 102 


104 



The sum of the Larmor and Thomas interactions takes the form 

- 4) ^i§^fW-WU^). (3.31) 

The factor g/m! — 1/m appears both in the spin-orbit and quadrupole-spin 
interactions. It gives the one-half reduction of the energy levels splitting for 
Dirac spin- 1/2 particles in heavy nuclei. 

Similarly, the interaction potential of the nucleus spin and the VL~ quadrupole 
moment can be found to be 

"■^ -miB-^) ^5(^-"(S.S)r-(I.L)fi;f.(0). (3.32, 

A modification of the contact interactions due to relativistic effects and fi- 
nite volume of the nuclei is discussed in Appendix C, with applications to the 
muonic atom. It is shown, in particular, that the L = 3 component of the elec- 
trostatic potential and the Darwin term generate contact P-wave interactions 
also. 
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3.5 Quadrupole-quadrupole contact P -wave interaction 



The contact part of the quadrupole-quadrupole interaction proposed in Ref. 
looks like 

Uq,q = -^Qz^Q^^lV^V^ - ^r^A)<5(x). (3.33) 

After integration over the angles and some additional algebra, one gets for 
L = 1 multiplet 



a 3Qz 3Q 
^QzQ 63 2/(2/- 1) 2S(2S- 1) ^ ' ^ 



X 



^r"^(I, I)r°^(S, S) - 3r^"(I, I)r"^(S, S)r^^(L, L) ) /^^^^(O). 



3.6 Numerical estimates of long-range and contact P-wave interactions in 
exotic atoms 



The contact P-wave interactions and the long-range interactions in ^^Nf2~ 
exotic atom are compared in Appendix B. In Tabled! we report the magnitudes 
of the interaction energies 

t/Pl = Uls + UjL + Uis + Uq,l + Ulq (3.35) 

and 

U^'^ = U^Q + Kzs + U^QzQ (3-36) 
for light nuclei with atomic numbers below 10 and spins / > 1/2. The quan- 
tities f/mal given in Table [1] for r = 1, 2 are defined by 



[2.] _ Tr[(^P^'] 
F \ Tr[l] 



(3.37) 



^mal is maximum over F-multiplets of the root mean square of eigenval- 
ues of f/P'']. f/J^[x t^lJlx typical strengths of the (f/c)^ and {v/cY 
interactions. We observe that f/^ax 2-5 orders of magnitude greater than 
[/W . 

^max 

One can see that the variations of < FJ'|f/W|FJ > in different nuclei are 
irregular and high. Counting of powers of f/c ~ aZ/n gives however the right 
first idea on the magnitudes of various interactions. 
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4 Decays of Q exotic atoms 



Decays of the Q~ exotic atoms proceed due to kaon t-cliannel exchanges be- 
tween nucleons of the nucleus and the ^"-hyperon: pQ^ A'EP + 180 MeV 
or ES + 100 MeV. These channels are shown on Fig. 2. 



Transition vertices: Relativistic expressions 



The effective vertices of the transitions p YK^ ,YK*^ where F = A, 11° 
and Q~K~^ Q^K*^ S*^ may be written in the form 

< Y\Jp{0)\N > = gNYKuiPY, SY)il5uiPN, sn), (4.1) 

< E\Jp{0)\n > = -^u{p^, s^)q^u^{pn, sn), (4.2) 

< Y\J^{0)\N >=U{PY, Sy){g^YK'Y + ^^^^cr^^'QuHpN, s^), (4.3) 

< E\Ji;{Om > = E ^^(PH, s^)r,^u,ipn, sn). (4.4) 

mi 

Here, Jp(x) is the pseudoscalar current coupled to the pseudoscalar mesons, 
Jy(x) is the vector current coupled to the vector mesons, and 



'^2yf, = -{,qvP,,- q- Pgvf,)!^, (4.6) 

The vertices of the transitions p S+iT", and n~K° S", n-K*° 

S~ are related by isotopic symmetry with fl4.ip - fl4.4l) . The amplitudes of the 
channels H^rP and are in the ratio —1 : ^/2. It is thus sufficient to 

calculate pQ' with Y = A, S°. 

The normalization and sign conventions of vertices fl4.ip and (14. 3 p follow to 
Refs. [nasi, the vertex (g^D is defined like in Ref. [Sj, the vertex (gaj) is 
simply related to that of Refs. [531168] . We use here dimensionless coupling 
constants, u and are relativistic spinors of the spin-1/2 particles and Q~ 
with the normalizations of Appendix D. Furthermore, 



q = Pn -Pe=Py -Pn, 

P = {pn+PN)/2 = ip3+PY)/2. {A.i 
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Fig. 2. Decay of pQ -exotic atom due to t-channel kaon exchange. 
4.2 Coupling constants 

The vertices —>■ 'EPK^ and pir^ are related by T±- and V±-spin 

operators of the SU3 symmetry group. The couphng constants for these chan- 
nels are equal in absolute values and opposite in signs. From the width of the 
A — > Nit decay one finds 

^ = -^ = 4.0. (4.9) 
v47r \/47r 

If the SU3 relations would include the mass parameter of the vertex (14.21) . the 
coupling constant Qq-bk could be m-^/mN = 1.4 times higher. It could result 
to an increase of the decay widths by a factor of two. 

Constraints from unitarity in the multichannel partial-wave analysis of the 
7iN scattering |66] are known to be not sufficient to find relative phases of the 
amplitudes with vrA^ and ttttN final states. The relative phase of Qan-k and 
dANp is thus unknown. 

The ratios of Qnak and Qqek, and Qnak* and fj^^lx* are known from the SUq 
symmetry (Appendix D). 

The coupling constants Qnak* and /nak* are related by SU3 to gNNp and 
Jnnp the relative sign of which is fixed by the interference of the vector meson 
exchanges in the nucleon-nucleon elastic scattering, by analytical continuation 
of the irN scattering amplitude to the t-channel [SZ] , and by the VMD model. 
The relative signs of f'^^p are fixed by the normalization to A photo- and 
electroproduction data [68]. The phase of the coupling constants does affect 
the interference between the K and ii'*-exchange diagrams. 

The vertex 'EPK*^ is related by SU^ with the vertex A++ PP^- 

The VMD model and isotopic symmetry allow to relate the latter with the 
A"*" p7 transition in which the magnetic form factor dominates. 

Like for pseudoscalar mesons, the vertices 'EPK*" and A++ pp^ con- 
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vert to each other by T±- and V±-spin operators. The corresponding couphng 
constants /^l/^. and /^Afp are equal in absolute values and opposite in signs. 
The vertex pp^ known from eVMD [HS] contains isospin factor ^2/3. 

In Ref. [68], the VMD model is extended in order to fulfil requirements of the 
quark counting rules by including higher radial states of the vector mesons. 
The ratios f^^Nv/dv are fixed for all radial states. The coupling constants 
/Aivv and gy are however known separately for the ground state p- and uj- 
mesons only. In our case, the transition momenta are not high, so we apply 
an integral description of the OBE amplitude by attributing the vector meson 
exchange potential to the ground state vector mesons. The most part of the 
experimental data come from the A electroproduction experiments i.e. from 
the spacelike region, so behavior of the transition form factors at = is 
determined more reliably than at = m?^. From static limit of the transition 
form factors, we obtain 

Jn^K' _ ^AA.p_2.73, (4.10) 
-1.68, (4.11) 
-1.42. (4.12) 



yivr \/47r 

f[2] /■[2] 

VAtt \/47r 

f[3] fm 

JHEK* _ JANp 



Air VAiT 

The pure magnetic transition would imply /'^l = — /'^l = — 2/^'^]. As shown in 
Appendix D, these coupling constants are real and their overal sign is fixed 
from the requirement that the p-meson coupling constants gNNp and /nnp are 
positive. If the ground state coupling constants are determined as residues of 
the transition form factors at = m^, one could get values 2.8 times higher. 
The final ambiguity should, probably, be smaller since radial states interfere 
at = destructively. In Eqs. (l4.10p - (14.12^ the SU^ symmetry is applied to 
the dimensionless coupling constants. If mass parameter of the vertex (14. 4p is 
included into the SU^ scheme relations, the coupling constants could increase 
by a factor of [ms/mN)^ ~ 2, while the contribution of the K* exchange to 
the decay rates could increase by a factor of four. 

The pseudoscalar coupling constants Qnak and Qnek are related by SU3 to the 
pseudoscalar coupling gNNw and to the F/D ratio. The i^-exchange contribu- 
tion to the pfl~ AS°, decays is determined by pseudoscalar couphng 
constants 

^^^^ -3.79 and ^ = 1.16. (4.13) 



/Arc VAir 
The quoted values are from Stocks and Rijken [51], model NCS97a. The K* 
exchange contributions to the pfl~ AS°, S°S° decays are determined by 
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vector coupling constants 



^ = -1.20 and ^ = -0.69, (4.14) 
l^ = -3.19 and = 0.32. (4.15) 

The model NCS97f of Ref. |M1 gives for f^Y-K* a value 3.5 times higher. 
This uncertainty does not affect the widths significantly, since the channel 
SS is not dominant. The uncertainties of other vector coupling constants and 
pseudoscalar coupling constants do not exceed 30%. 



4-.3 Transition vertices: Non-relativistic reduction 



We describe decays to the final state A'EP. The decay channel is distinct 
by coupling constants and masses of the involved particles only. The decay 
amplitude for the channel S+S^ is known from the isotopic symmetry. In 
what follows, we work in the rest frame of pQ" . 

The energy released in the decays is small, so we apply nonrelativistic approx- 
imation. The Coulomb interactions of the charged virtual kaons are neglected. 
The nonrelativistic reduction of vertices (14.11) - (14.41) gives 



gNAKu{pA, SA)i-f5u{pN, Sn) = (pXCiQ^ipN, (4.16) 

u{Ps, s^)q^Uf,{pn, Sq) = '^tC2q'Cn^ (4.17) 



m(Pa, SK){,gNKK*l^ + ^^^ia°''q„)u{pN, Sjv) = v^aCsv^jv, (4.18) 

2m AT 

«(PA, SA)i9NAK'l°' + ^^^'>'(^°"'qu)uiPN, Sn) (4.19) 

2mN 



^i(C4pS + Cfg^)^ 



E ^^w(pH, s^)Tlu,{pn, sn) (4.20) 



E ^Hpb, ss)Ttu^{pn, sn) = ^tCf^^Cn, (4.21) 



where 
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= 


2mN ' 






Cs = 


Onak* , 




9nak* 


rriN 




7^ {-gNKK*^^ 


Og — 




^afi-y 


■Jqek' "t Jqek* 




— I ;^ 1 

2m| 


Og — 





(4.22) 

(4.23) 
(4.24) 
(4.25) 

+ {9NKK* + W*)^e°^^^^), (4.26) 
(T^(5°^), (4.27) 

ct"5^^, (4.28) 

(4.29) 

Here, the terms (m^ — m-^)/m-^ ~ (mA — 'mN)/'mN ~ q^/m|^ ~ ("^/c)^ are 
systematically neglected. 

In general, the iVA vector current is not conserved. Its longitudinal component, 
however, does not contribute to the decay amplitude, since the divergence of 
the flS vector current defined by (14. 4p vanishes both on- and off-shell. It 
allows to calculate the K* exchange amplitude as the product of A^A and fiS 
vector currents. In the limit of = the current (14.31) is conserved, so 
its divergence ~ (mA — ■mi\i)/m]\; can be neglected in the nonrelativistic limit 
anyway. 



4-4 Decays of pVt atoms 



The pVt atomic state is defined by 



\P,SN3,sn3,nlm >= J ^^^^3 d{P-pN-pn)^nim{ )|Pjv, sa^s; Pn, sm > 

The final state consists of a plane wave of A and S: 

|pa,sa3; Ph, sh3 > • (4.30) 



These states are normalized by 
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< P',s'jv3, Sm, n'l'm' I P,SAr3, sqs, ulm > 



= (2vr)35(p;-pJ(27r)35(p^ - Ps)Ss',^,s,.A^„s-^s- 

The first condition is a consequence of tlie ortliogonality of liydrogen-like wave 
functions 

^^;^*n'Z'm'(q)*^nim(q) = Sn'Jl'lSm'm- (4.31) 

Consider the lowest order S"- matrix element responsible for decay pQ" AS*^: 



Sfi = f J J d^xd% < pA,SA3;p~,Ss3\T[CEQK{x)CANK{y) + 

Entering this expression 

Cank{x) =(f)^{x)JAN{x), (4.32) 

C^nKix) = Jsnix)(f){x), (4.33) 

/:ANK*{x) = -<Pti^)JM, (4.34) 

^snK*ix) = -J^^{x)(l)f,{x). (4.35) 

are effective Lagrangian densities corresponding to vertices (14. ip - (14.41) . 0(a;) 
and (pfxix) are K~^- and i^'"^*-meson fields. 

In the nonrelativistic approximation, the amplitude takes the form 



^f^=J:'^f^ (4.36) 

where 



A 



m], = -HCMHC^J^'Cnl (4.37) 
^% = [vICsVn] MiCfVsJy + Cf'jf)enl (4.38) 



= -H{C,P^4 + Cf4')^NMCt'Cn]- (4.39) 



The first matrix element corresponds to the /T-meson exchange, the last two 
ones correspond to the i^'*-meson exchanges as the products of timelike and 
spacelike components of the transition vector currents. 

The functions Jm "'^" entering Eqs. (14.391) s = 1, 2 appear upon the integration 
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over the atomic wave function: 



-(PH) = /*^5(p.+p„)*„45^)^,« (4.40) 

where q = pa — Ps = Pa — Pn and M = P,V. These functions are symmetric 
in ai...as and can be converted to the form 



/_• _ d d e~"^*f^^ 



Here, use is made of the momentum conservation Pat + = Pa + Ps = 0. 
The bound state wave function is written then in the coordinate representa- 
tion. The effective meson mass entering the Yukawa potential equals = 



^/iTT'M - Qo where 



2 V - V " 2 "7 ' " 2v^ 



and a/s = mn + rriN + En- The effective masses of K- and ii'*-mesons, m*i^ = 
410 MeV and m^, = 850 MeV, are distinct from the vacuum masses rriK — 494 
MeV and mx* = 892 MeV. 



We wish to find 



Y: \mM' = YnAB, (4.41) 

spins A,B 



where 



spins 



The straightforward calculations give 
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ni2^n;, = o, (4.43) 

—n^*— ^ 9nak gnEK Qnak* + Inkk* T^ufn^K* 

'MS — — 2 

X ( J°V^^* - Jr^y^*), (4.44) 

o 4 ~ 

6 m= 



+ + /ili.0'(3J^'J^'* - Jr^y'*)), (4.45) 

^23 = 7^;2 = ^"^f (p| J^- - pg J^pI J^-), (4.46) 

— 'Z^y^9N\K*\PB.'Jv'Jv - Pe'^vPe'^v ) 

HgNKK* + fNKK^f{J^^jf* + J^Jv*))- (4.47) 



The values T^n, 7^22, "^23, and 7^33 are not affected by phase ambiguities of 
the couphng constants. The phase factor of the interference term TZ^^ is fixed 
in Appendix D. 



The functions J^'""' for s = 0, 1, 2 can be expressed as follows: 



JMip)^i-iyYU^)jSip), (4.48) 
P 

JM{p) = ^i-iyyim{^)jS{p), (4.49) 

+\s'^^{-iyYU-)J^S{p) (4.50) 

where 
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jSiP) = - / r^drRr,i{r)ji{pr) , (4.51) 

= - / r'drR^r) (Mpr) - 5,o) + - , (4-52) 

JO y J, ^^J J. 

Jm{p) = " r'^drRniir) [ji{pr) - (1 - ^pV^)5io - ^J5r(5;i^ 

xf!f (4,53) 

\ ,y-iO / 

= -<iJtKp) - RnliO). (4.54) 



It is customary to remove 5(x) terms originating from meson exchanges from 
baryon-baryon potentials, since such terms are removed anyway by zeros in 
the baryon wave functions in presence of the repulsive core. We thus drop the 
Rni{0) term in Eq. fH^Ijl . 

To arrive at Eqs. (14.511) - (I4.54p . we used decomposition of the plane wave 

e-'^- = A7rY.^-'Mpr)YU-)Y;J-) 
t: p r 

and decomposition of the atomic wave function 

^n«m(x) = Rnl{r)Yirn{-). 



The integrand of j\J (p) is regularized to insure the possibility of interchanging 
order of the integration over r and the differentiation over p". The replace- 
ment ji{pr) — > ji{pr) — 6io does not affect the result, since for / = J^m{p)' 
enters Eqs. (l4.55|) - (I4.58P only. The integrand of (p) is regularized at r = 
similarly. The functions J^^ (p) and j|f (p) are well defined. 

The functions (I4.48P - (14.501) have dimensions (in units h = c = 1): [Jm'""" (p)] = 
[M-V2+^], and also [jfjip)] = [M-'/% [jg(p)] = [M^^ [jg(p)] = [M^^ 
and [jS(p)] = [M3/2]. 

By integrating the decay amplitude squared over ps, four integrals appear 



4L'(P) = / dn^J2hp)J^^,*ip), (4.55) 

41m'(p) = / dn^J^-{p)J^Yip), (4.56) 

KfkAp) = J dn^p'-JMp^Jlm. (4.57) 

k^LAp) =p'I dnMP)JM'iP)- (4.58) 



22 



Table 2 

Decay widths of atoms [in keV] in the lowest atomic levels n = 1 4 and 
L = 0-^3. The results are given for the AH and decay channels. The last 
channel includes the summation over the TPrP and states. 



L 





1 


2 


3 


n channel 


r X 10^ 


r X 10** 


r X 10^^ 


r X lo^^ 


1 AH 


60 








SH 


1.8 








2 AS 


7.5 


4.1 






ss 


0.0 


2.2 






3 AS 


2.3 


1.4 


2.0 




ss 


0.0 


0.8 


1.7 




4 AS 


1.0 


0.6 


1.2 


1.0 


SS 


0.0 


0.4 


1.0 


0.6 



In Appendix E, we provide identities useful for the integration over the angles. 
The integrals fl4.55p - (14.581) are found to be 



-iA,j£l(p)A,jg,(p) + ijg(p)jSl(p), 
K^SM'ip) = J^kp)jS'ip), 

K^SM'ip)=p'Jikpyjikpy, 

Kf,M'{p) =p'j^Kpyjlkpy + i{i + i)4kp)4/kp), 

where 



The atomic decay width averaged over the spin states of the proton and VL~ 
equals 

where 



fdQ^^Y.'^AB, (4.59) 

A,B 



J{s- (mA + m^f) (s - (mA - 7713)2) 
PS = (4.60) 

is the center-of-mass H-hyperon momentum. In Table [2] we give decay widths 
of lowest pVt'^ atomic levels. 

The total decay width to the channels SS is enhanced by isotopic factor 3. In 
P states, the interference of the K and K* exchnage amplitudes is destructive 
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for the AS channel and constructive for the ES channels, producing thereby 
an additional enhancement of the SS channels. For 2P state, we get Fah = 
(6.1 + 2.3 - 4.4) X 10-6 keV and Fsh = (1.8 + 0.2 + 0.3) x 10"^ keV. The first, 
second and third numbers correspond to K, K* and their interference. The 
interference pattern of the 5- wave is opposite. It is discussed in Sect. 4.5 and 
Appendix D. 

Since Vl^ spin is high, one can expect nontrivial F-dependence of the decay 
widths. The S-wave decay from an F = 2 spin state e.g. can be suppressed, 
since it is coupled to a D-wave channel only. 

In Eqs. fl4.5ip - fl4.54p one can expand Rniir) around r = and keep the lowest 
order term due to the small binding energy. As a result, the widths scale with 
n approximately as squares of the L-th order derivatives of Rnbif) at r = 0. 

The decay width of the 2P state is one order of magnitude lower as compared 
to pp atoms |i55j. The pp atoms are coupled to the continuum through the 
annihilation channel only. p'E~ and pS" atomic states decay due to t-channel 
kaon and pion exchanges [TUfSS] . Such decays can be calculated like decays of 
pVL~ atoms. 

The formalism presented here is similar to that used in Ref. [5S] for calculation 
of weak decays of a loosely bound hypothetical H particle [57] , which proceed 
due to kaon exchange between A-hyperons. 



4-5 Decays of high-Z nucleus-Vt atoms 

Strong decays of Vt~ exotic atoms with high-Z nuclei proceed under kinemat- 
ics conditions which are more complicated as compared to pVL~ decays. A 
microscopic approach to calculate the strong decays of hyper on atoms is 
discussed by Loiseau and Wycech [58] . It is based on the impulse approxima- 
tion where it is assumed that in the final state one has a nucleus left with a 
hole in a single particle state and two particles as plane waves. 

The optical potential method represents a conventional approach to calculate 
strong decay widths of high-Z exotic atoms. It consists in the determination 
of free scattering lengths and finding the average value of imaginary part of 
the nuclear optical potential. Such a method is in the qualitative agreement 
with the decay rates of S~ exotic atoms [59] . 





The imaginary part of the NVt scattering lengths 



aj = lim 



(4.61) 



4:7i{ma + miy) 
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can be estimated using the plane wave Born approximation (PWBA) from the 
OBE diagram of Fig. 2: 

ai = 2.2 fm (4.62) 

both for proton and neutron. This value has the content a/ = 0.45 + 1.31 + 
0.09 + 0.03 + 0.44. The first and third numbers come from the i^'-exchange, 
the second and fourth numbers come from the if*-exchange, the first two 
numbers and the next two ones correspond to AH and SS channels, and the 
last number is the interference. 

As proposed by M. Ericson and T. E. O. Ericson [60], the atomic decay widths 
can be calculated in terms of the scattering length aj and integral over the 
nuclear density 

47]- nR 

r = —ai / r^drp{r)RlL{r), (4.63) 
/i Jo 

where R = ^ < >^/^= r^A^^^ is the nucleus radius. For 12 < A < 40, 
To = 1.35 fm and for heavy nuclei tq = 1.20 fm [TD]. The density p(r) is 
normalized according to / p{r)dV = A. The value of is the reduced mass of 
Q~ and nucleus. 

The estimate (14.631) neglects modification of the wave function due to finite 
nuclear volume, relativistic effects, and real part of the strong interaction 
potential. In order to remove some of systematic errors caused by neglecting 
those effects, we normalize the decay width to F = 0.05 keV of ^^NS~ exotic 
atom in 3-D state calculated by Batty, Friedman and Gal [15]. Using Eq. (l4.63p . 
we reproduce then with reasonable precision the reported theoretical values 
of the decay widths of S~-atoms formed with ^^C, ^^O, and ^^F in 3D states 
and ^^Si in 4F state. 

The pseudoscalar coupling constant Qsak/V^ = 1-10 [54j|, model NCS97a, 
is small as compared to gnsx/V^ = —4.0 given by Eg. (14.91) . The channel 
SA where the pseudoscalar coupling constant gssx/V^ = —4.69 is large 
is blocked by the energy conservation. One can expect that the K exchange 
contribution to the pQ~ scatering length is higher as compared to that of the 
A^H~ scattering length. The amplitudes pQ~ — > AS due to the K and K* 
exchanges are in the ratio 1 : 1.7. 

The imaginary part of the NE~ A A scattering length is only a/ = 0.04 fm 
|15j . This small value is the result of small coupling constant q-ekki statistical 
weight 1/16 of the AA channel, and small phase space of AA. Increasing the 
imaginary potential causes the calculated widths to increase by roughly the 
same proportion. Taking into account the decreased Bohr radius due to the 
higher mass of VL~ and significantly increased imaginary part of the scattering 
length, we obtain decay widths of the Vt~ exotic atoms reported in Table [1] 
Using the same assumptions, we obtain for ^°^Pbfi~ atom Fg ~ 7 keV in 
L = n — 1 = 9 state, Fio ~ 0.1 keV in L = n — 1 = 10 state, and Fn ~ 2 eV 
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for L = n — 1 = 11. In Table [H we restricted by giving estimates for nuclei 
A > 6, since Eq. fl4.63p matches smoothly at A = 1 with decay widths of S 
levels, whereas we list decay widths of P levels. 

The G-matrix formalism by Yamamoto et al. |61| results to decay widths of H 
atoms one order of the magnitude higher [15]. Accordingly, one can consider 
above estimates of widths as the lower bounds. 

The value of the scattering length Eq. fl4.62l) is three times higher as compared 
to the pK~ scattering length. Uncertainties in the coupling constants of 
discussed in Sect. 4.2, point in many cases towards even higher values of ai. 

We see from Table [T] that widths of the 2P states are several orders of the 
magnitude higher than those required to differentiate the contact P-wave in- 
teractions and up to three orders of magnitude higher than the long-range 
interactions. The kaon exchange in circular orbits of the ^°^Pbi7~ exotic atoms 
has negligible effect starting from L ~ 10. 

The strong-interaction shift is usually expected to be as large as the width 
and a fraction of the strong-interaction shift could also be spin-dependent, thus 
contributing to the apparent hyperfine splitting. Lorentz scalar component of 
the two-kaon exchange potential does not generate hyperfine splitting. The 
vector component generates hyperfine splitting due to spin-tensor interaction 
AE ~ {otZY X 1 MeV ~ 3 keV (for mass and spin dependence of OBE 
potentials of spin-1/2 particles see Ref. [13], the numerical values are quoted 
for ^^NVt^). The effect is comparable with the long-range quadrupole-orbit 
interaction. 

In Ref. [7], the circular transitions in ^°^Pbr2~ exotic atoms for L = 10 — > 9 
levels are estimated. The transition energy is about 0.5 MeV, each level splits 
into four sublevels due to spin-orbit (~ 2 keV for L = 9) and quadrupole-orbit 
(~ 0.2 keV for L = 9 and Q ~ 3 x 10~^ fm^) interactions. In the experiments 
with atoms [2], the peak positions of the photon energies are determined 
with accuracy of a few tens of eV. Using the same technique, it would be 
possible to measure the fl~ quadrupole moment with an accuracy of ~ 30% in 
the circular transitions L = 11 ^ 10. The strong interaction shift of the lower 
L = 10 level is expected to be ~ Fio, while splitting due to strong interactions, 
sensitive to the quadrupole moment, is smaller: AE ~ {aZ/nyViQ < 1 eV. 

In Appendix C, we compare the natural widths and the contact P-wave inter- 
actions in /i-meson exotic atoms. The estimates reported in Table IC.ll suggest 
that in high-Z nuclei the contact interactions are of order of the radiative 
widths of 2P3/2 levels. 
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5 Conclusions 



In this work, we investigated the possibihty of measurement of the quadrupole 
moment by observing X-rays from low-L transitions in Q~ exotic atoms. The 
magnitude of fine and hyperfine sphttings of the energy levels has been com- 
pared to strong decay widths caused by reactions NQ~ — >■ AS, SS in Q~ atoms 
formed with light stable nuclei with atomic numbers below 10 and spins above 
1/2. 

We proposed, firstly, a minor modification of the Q~ spin-orbit interaction 
used in the earlier works [T][T7] in order to bring it in agreement with theory 
of Thomas precession. 

Secondly, we described new kinds of the contact P-wave interactions for par- 
ticles with electric quadrupole and magnetic dipole moments. We found that 
Thomas correction for the quadrupole- spin contact P-wave interaction is the 
same as for the spin-orbit long-range interaction. The Darwin term connected 
to the particles Zitterbewegung represents yet another source of contact P- 
wave interactions. 

The long-range interactions appearing to the order (f/c)^ such as spin-orbit 
interactions, spin-spin tensor interaction, and quadrupole - orbit interactions, 
have been discussed and included into the numerical estimates of the energy 
splitting. We showed the for 2P states of Q~ exotic atoms with stable nuclei 
up to ^^F the contact interactions are 2 — 5 orders of magnitude weaker than 
the conventional long-range interactions. The quantitative evaluation of the 
contact P-wave interactions suffers from the poor knowledge of the short-range 
component of the baryon-baryon interactions. 

The contact P-wave electric quadrupole - magnetic dipole interaction exists 
in ordinary atoms and /x-meson atoms. In high-Z nuclei, the magnitudes of 
the contact P-wave interactions in /x-meson atoms are comparable with the 
natural widths of 2P3/2 levels. 

Thirdly, we calculated strong decay rates of pfl~ exotic atoms due to reactions 
pQ^ AH, SS caused by K and K* t-channel exchanges. The problem is 
solved analytically for arbitrary principal and orbital quantum numbers. The 
decay rates of the lowest pQ~ atomic levels, averaged over the proton and Q~ 
spin states, are reported in Table [2l 

Rough estimates of strong decay rates of the fl~ exotic atoms formed with the 
light nuclei and of the ^"^Pbfi" exotic atom have been made. For 2P states of 
the \ow-Z Q~ atoms, we get strong decay widths up to three orders of magni- 
tude higher than splitting caused by the conventional long-range interactions 
and, respectively, 5-6 orders of magnitude higher than the contact P-wave in- 
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teractions. Table [T] summarizes the estimates of magnitudes of the interactions 
and the strong decay widths. The contact P-wave interactions are not useful 
for measurement of the Q~ quadrupole moment. 

Strong decay channels in ^°®Pbfi" exotic atoms are small in the circular tran- 
sitions starting from L = n — 1 ~ 10. Such transitions minimize theoretical 
uncertainties inherent to the problem and can be suitable for measurements 
of the fl~ quadrupole moment. 

The authors are grateful to M. D. Semon for correspondence and B. V. Marte- 
myanov for discussions on the Thomas precession effect. The authors wish to 
acknowledge referees of Nuclear Physics A for useful remarks and suggestions. 
This work is supported by RFBR grant No. 06-02-04004 and DFG grant No. 
436 RUS 113/721/0-2. 



A Thomas precession 



A formally complete treatment of the contact P-wave interactions would re- 
quire the knowledge of other interactions and, specifically, the rate of Thomas 
precession to order (f/c)^. There exists a controversy in the evaluation of 
Thomas precession effect beyond iy/cY as discussed recently [72]. Here, we 
give relativistic treatment of Thomas precession. 

Let us consider coordinate systems K and K' in which a particle has four- 
velocities u = (7,7v) and u' = (1,0), respectively. In the coordinate system 
K, particle moves with velocity v, whereas in the coordinate system K' it is 
at rest. Given that w = {wo,w) in K, we search for the coordinates of w' in 
K'. 

We split w into the parallel and transverse components with respect to the 
direction of the velocity: W|| = n(nw) and = w — n(nw), where n = v/f . 
The transformation of {w^,w\\) where w\\ = (nw), are well known while w_|_ 
does not transform. One can write therefore 



w 



WQ = 'y{wo - vnw), 
nw' = 7(nw — vwq), 
nfnw') =w — nfnw). 



(A.l) 
(A.2) 
(A.3) 



This system of equations allows to find the Lorentz transformation matrix (see 
e.g. |Mj) 



7 V7 

V7 1 -|- (7 — l)n 



(A.4) 



n 
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The coordinates of a four-vector w in two Lorentz coordinate systems K and 
K' are related by 

L(v)w;' = w. (A.5) 
Remind that K' moves with velocity v in K. 

The particle polarization is a three-dimensional unit vector, a. It is defined 
in the particle rest frame. Relativistically, polarization is characterized by a 
four-dimensional vector a. Given in K a particle with a four-velocity u and a 
polarization four-vector a, such that 



m2 = 1, (A.6) 
a^ = -l, (A.7) 
a-u = 0, (A.8) 

one can define a three-dimensional unit vector, a, as space-like component of 
a' from equation 

L(v)a' = a, (A.9) 



where v is the particle velocity in K. We thus make boost to K' where 
u' = (1,0), a' = (0,a). Three-dimensional vector a is called polarization of a 
particle which moves with velocity v in K. 

Let us consider particles 1 and 2 with four- velocities ui and U2 and polarization 
four-vector ai and 02 in K. They can be considered as two different states of 
the same particle, separated by a time interval 6t in K. Particles 1 and 2 are at 
rest in coordinate systems K' and K", respectively, as shown on Fig. lA.ll As 
discussed above, we can define two three-dimensional unit polarization vectors 
ai and a2: 



a; = (0,ai) =L(-vi)ai, (A.IO) 
a'2' = (0,a2) =L(-V2)a2, (A.ll) 

velocities Vi and V2 are known since Ui and U2 are known. 

Four-velocities U2 and Ui are related by a Lorentz boost. We denote K'" a 
coordinate systems obtained from K' by such Lorentz boost: L{5\)w"' = w'. 
In particular. 



L(5vX' = m;, (A. 12) 

L(5vX = u'^. (A. 13) 

Particle 2 is at rest both in K" and K'", K" and K'" are related by a rotation. 
Taking into account that u'2 = u[ = (1,0), we get 

L{6y)u[ = u'^. (A.14) 
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K",K"' 




Fig. A.l. Schematic representation of Lorentz boosts relating coordinate systems 
K, K', K" , and K'" involved into calculation of the precession rate of polarization 
of a particle moving along a trajectory. 

We require that polarization four-vectors be related by the same transforma- 
tion: 

L(5v)a; = a2. (A. 15) 

This ensures to fulfil Eqs.(jAj]) - flAlSl) for particle 2 provided Eqs.dXl]) - flAlSl) 
are fulfilled for particle 1. The way the four- vectors are related defines parallel 
transport from K' to K'": The coordinates of all four- vectors attributed in K' 
and K'", respectively, to particles 1 and 2 remain unchanged. In particular, 
a'2" = a; = (0,ai). 

The relativistic composition of velocities can be used to express V2 in terms 
of Vi and 5v: 

V2 = vi © (5v = vi + (- — r^i ® vi)(5v, (A.16) 

7 7 + 1 

where 7 = 1/ yl — vf . According to an observer in K'", K moves with velocity 
V3 7^ -V2: 

- V3 = 5v© vi = vi + (1 - vi (g) vi)5v. (A.17) 

The composition of velocities is defined by the composition of Lorentz boosts: 
U2 = L{yi)L{5w)u2 = L{-Vi © 6\)u'2. One finds that V2 = vi © 5v is velocity 
of K'" in K. K" and K'" are distinct by a rotation, so V2 is velocity of K" in 
K too. 

Now, it is straightforward to find 



a'2 = L(-V2)a2 
= L(-V2)L(vi)a2 
= L(-V2)L(vi)L(5v)a'2" 

= L(-V2)L(vi)L((5v)a;. (A.18) 
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These equations show that 

a2 = Mai, 

where M is a rotation matrix, such that 



1 



L(-V2)L(vi)L((5v). 



(A.19) 



(A.20) 



Applying R, one gets 
where 



(5ai 



Lai = ai + SsLi 

7 



7 + 1 



vi X 5v) X ai 



(A.21) 
(A.22) 



Circular motion where Vi^v = implies 

Sv 



7 



CJ X Vi5t, 



(A.23) 



where cj is the orbital rotation frequency, 6t is a time interval in K, and 
(5v/7 = V2 — vi according to Eq. (lA.16p . so one gets 



6SLi = (1 — 7)cl» X SLi5t. 



(A.24) 



This equation shows that vector ai experiences a precession in K with fre- 
quency fix given by Eq. fl3.10l) . 



In Refs. 



Qt is 7 times smaller. The possible reason of the discrep- 



ancy might be the noncommutativity of relativistic composition of velocities. 
The velocity of K'" in K and minus velocity of K in K'" deviate from Vi for 
circular motion by 6\/'~f and 5v, respectively. The quantity a; x \i6t refers to 
the variation of velocity of K'" in K. 



The definition of polarization of a moving particle using Eqs. (1A.10|) and (lA.lip 
allows to attribute the transparent physical meaning to Lorentz boosts relating 
the coordinate systems K . . . K'" . In all coordinate systems related by Lorentz 
boosts with the particle rest frame we observe the same three-dimensional unit 
polarization vector. If, however, a particle has been accelerated by a sequence 
of non-coUinear Lorentz boosts, its polarization does rotate. The particle polar- 
ization and, accordingly, its time evolution depend on the coordinate system. 
From the point of view of an observer in K' there is no rotation in transit 
from K' to K'" . However, in K we do observe a rotation. 

Relativistic expression for Thomas precession frequency an external electro- 
magnetic field can be found using Eq.(A.22) with vi x 5v = 7V1 x (v2 — vi) = 
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7Vi X \iSt and the covariant equation of motion for charged particles: 



= vi X (E + vi X B). (A.25) 

7 + 1 m 



All quantities entering this equation are defined in the laboratory frame K. 
The Larmor precession frequency in the co-moving frame K' can be found 
using the Lorentz transformation for electromagnetic field: 



n'^ = -|b' = -^(7(B - vi X E) - (7 - l)ni(ni • B)). (A.26) 



Here, rii = vi/|vi| and B' is the magnetic field in K' (cf. Eq.(III.8)). The sum 
fix + ^l/i gives the total spin precession frequency in the laboratory frame, 
in the exact agreement with the Bargmann-Michel-Telegdi equation [H] . The 
above arguments do not rely on the assumption of = 1/2. 

Equation flA.24p is in agreement with Refs. [7Hl75|76] . It is physically equiva- 
lent to the equation for rotation of axes of a "Born- rigid electron" on circular 
orbit, derived first by Foppl and Daniell [77]. 



B Matrix elements of angular momentum operators and fine and 
hyperfine splitting in ^^Nfi" exotic atom 



Let us consider a symmetric tensor r"'^(a, b) constructed in terms of operators 
a and b: 

r"^(a, b) = a°6^ + a^6" - -a • b(5"^. (B.l) 

3 

In our case, a, b = F, I, J, L, S where F = I + J is the total angular 
momentum of the system and J = L + S is the total angular momentum of 
n~. Recall that [a",a'3] = ie^^^a^ for a = F, I, J, L, S, [F",a^] = ie"^^a^ 
for a = I, J, L, S, [J", a^] = ie^^^'a^ for a = L, S, in other cases [a", I/] = 
for a 7^ b. 

Let a + b = c. Consider contractions of two tensors r"^: 
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z/(a,b,c)=r°^(a,b)r"^(c,c 
ta ■ cb ■ c — 
z/(a,b)=r"'^(a, a)r"'3(b,b) 
la ■ ba ■ b + 

z/(a,c)=r°^(a, a)r"'^(c,c) 



4 

4a ■ cb ■ c - -a • be • c, (B.2) 



4 

4a ba b + 2a b a ■ ab • b, (B.3) 

3 



4 

= 4a ■ ca ■ c — 2a • c — -a ■ ac ■ c, (B.4) 
z/(c)=r"'3(c,c)r"^(c,c) 

= ^c(c+l)(2c- l)(2c + 3), (B.5) 



where c • c = c(c + 1). 
Specifically, we define 



[n, L) = r"^(n, n)r"^(L, L) = .l^^^^. (B.6) 



Upon averaging over states with fixed L (first line) or J (other lines) one can 
write 



r"^(n,n) 



KL) 


(B.7) 


z/(J) 


(B.8) 


KS,L,J) 
KJ) ^ 


(B.9) 


z/ J 


(B.IO) 



r"^(S,S) 
r"^(S,L) 
r"'3(L,L) 



The contraction of three functions r"^ entering Eg. (13.341) . averaged over a 
fixed J state, gives 



/.(S, L, J) = r"'3(S, S)r^^(L, L)r^"(J, J) 

= (2L ■ S - i)z/(L, S, J) - h{S + l)z/(L, J) - ^L(L + l)z/(S, J). (B.ll) 



Using Ref.jH], one gets 



1 _ 2 
^ ^ ^~ n3L(L + l)(2L + l)a| ^ 
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and, for P-wave, 



5 ' 
B 



where as = l/{aZm') is the Bohr radius. 

The diagonal matrix elements of the interaction energies in the J' J basis can 
be found to be 



- "4 2/(2/-l) < ^ KlRJ^ ' ^ ^ 

3Q ^ 1 ^ z/(n,L)z/(L,S) 
= -~2S{2S-1) < ^ > ' ^^-^^^ 



rrcL 3Qz t/(I,J)t/(S,L,J) 

^^-^ - 40mm' 2/(2/ -1) ^^"^^^ 



U' 



cL 



aZgz SQ I ■ J 



~ lOMm' 25(25- 1) J- J 

a 3Qz 3Q z/(I,J) 



X (^L-SJ-S- ^L- JS-S)/?;"i(0), (B.18) 



63 2/(2/- 1) 25(25- 1) z/(J) 

X (^LiS, J) - 3^(S, L, J)) R'^M- (B.19) 

Other matrix elements can be calculated using elementary tools. 

The diagonal matrix element of Uq^^l in the basis of fixed I + L has the 
form of Eq. flRTSj) with z/(L, J)i/(I, J)/z/(J) replaced by i/(I, L) (cf. Ref. [M], 
Chap. VI). In J = 1/2 states of the Q~ atoms, the diagonal element of the 
quadrupole - quadrupole interaction (1B.19P vanishes, since J = 1/2 states 
do not have quadrupole moments. The diagonal matrix element of Uq^q for 
J 7^ 1/2 is calculated in Refs. p^lITT] . Equation (IB.19P is in agreement with 
Ref. 117]. 

We use for calculations of the Clebsch-Gordan coefficients a code provided by 
Sierra ^63j. The diagonal matrix elements for the potentials entering ^/'^l and 
t/'^l are in agreement with those calculated numerically. 
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Table B.l 

J' J matrix elements of fine and hyperfine interactions of order {v/c)"^ in the ^^Nfi" 
exotic atom for the n = 2, L = 1 state. LS stands for the interaction energy Uls 
Eq. (j3.13p . IL stands for the interaction energy Ujl Eq. (|3.14p . and so on. Parameters 
used in the calculation: nz = 0.404 n.m. [46j, ^ = -2.02 n.m. [311], Qz = 2.00 fm^ 
\0,Q = -2.8 X 10-2 fm^ [23]. The energy is given in keV. 



UM, keV 


LS 


ILxlO 


ISxlO 


QzL 


LQ 


F 


J 


1/2 3/2 5/2 


1/2 3/2 5/2 


1/2 3/2 5/2 


1/2 3/2 5/2 


1/2 3/2 5/2 


1/2 


1/2 
3/2 


-0.70 0.00 
0.00 -0.28 


0.11 0.12 
0.12 -0.11 


0.14 0.20 
0.20 -0.46 


0.00 1.51 
1.51 -2.70 


-0.33 0.00 
0.00 0.26 


3/2 


1/2 
3/2 
5/2 


-0.70 0.00 0.00 
0.00 -0.28 0.00 
0.00 0.00 0.42 


-0.05 0.20 0.00 
0.20 -0.04 0.12 
0.00 0.12 -0.23 


-0.07 0.31 0.00 
0.31 -0.18 -0.16 
0.00 -0.16 0.17 


0.00 -0.48 3.51 
-0.48 2.16 -3.48 
3.51 -3.48 1.89 


-0.33 0.00 0.00 
0.00 0.26 0.00 
0.00 0.00 -0.07 


5/2 


3/2 
5/2 


-0.28 0.00 
0.00 0.42 


0.07 0.15 
0.15 -0.65 


0.28 -0.20 
-0.20 0.05 


-0.54 1.86 
1.86 -2.16 


0.26 0.00 
0.00 -0.07 


7/2 


5/2 


0.42 


0.16 


-0.12 


0.68 


-0.07 



The contact interactions contribute to splitting of L = 1 states and mixing 
and splitting of L = and L = 2 states. 



The numerical magnitudes of the contact P-wave interactions and the long- 
range interactions are compared by considering splitting of 2P energy levels 
of the ^^Nfi- exotic atom. The nucleus ^^N has spin 1 = 1 and, respectively, 
magnetic and quadrupole moments. 

The matrix elements of fine and hyperfine interactions are calculated in the ba- 
sis J' J at fixed F for eight different contributions: spin-orbit LS and IL, spin- 
spin IS, quadrupole-orbit QzL and LQ, quadrupole-spin QzS and IQ, and 
quadrupole-quadrupole QzQ interactions. The results for various terms enter- 
ing the potential fl3.35p are shown in Table [RT] and for the interaction (13.361) in 
Table IB.2I We restricted ourselves with estimates of the Larmor components 
of the quadrupole-spin interactions. In 2P state, the matrix elements of the 
contact interactions < FJ'|f/W \FJ > are suppressed as {aZy ~ 3 x lO^'^ with 
respect to the matrix elements of the long-range interactions < FJ'\U^'^'^ \FJ >. 

Our estimate of the contact P-wave quadrupole-quadrupole splitting in ^^Nfi" 
is two orders of magnitude smaller than the estimate reported in Ref. [17] . 
The charge radius included into the estimate of Ref. [IT] can increase the 
hyperfine splitting, since Qq has the smallness ~ v/c or even ~ [v/c]^ as 
compared to the proton and charge radii, as discussed in Sect. 2. In Ref. 
[T7] . it is assumed that has the same magnitude as Qn, so the reason for 
the discrepancy is unclear. 

^ This manuscript was in press when the authors were informed that G. Karl and 
V. A. Novikov revised their estimate of the hyperfine splitting. Their new estimate 
is in agreement with our [V. A. Novikov, private communication]. 
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Table B.2 

J' J matrix elements of contact P-wave interactions of order {v/c)^ in the ^^Nfi" 
exotic atom in the n = 2, L = 1 state. IQ stands for the interaction energy Ujq 
Eq. (j3.32p and so on. Parameters and notations are the same as in Table [BTTl 



UW, keV 


IQxlO* 


QzSxl02 


QzQxlO^ 


F 


J 


1/2 3/2 5/2 


1/2 3/2 5/2 


1/2 3/2 5/2 


1/2 


1/2 
3/2 


-0.26 -0.03 
-0.03 -0.21 


0.00 -0.73 
-0.73 0.75 


0.00 0.10 
0.10 -1.02 


3/2 


1/2 
3/2 
5/2 


0.13 -0.05 0.00 
-0.05 -0.08 0.09 
0.00 0.09 0.11 


0.00 0.23 -0.49 
0.23 -0.60 -0.24 
-0.49 -0.24 0.78 


0.00 -0.03 -0.96 
-0.03 0.81 -0.37 
-0.96 -0.37 -0.07 


5/2 


3/2 
5/2 


0.12 0.11 
0.11 0.03 


0.15 0.13 
0.13 -0.90 


-0.20 0.20 
0.20 0.08 


7/2 


5/2 


-0.08 


0.28 


-0.03 



The spin-orbit interaction is not dominant, probably except for F = 1/2, 
so the total Q~ angular momentum J does not provide diagonal basis. For 
F = 3/2, 5/3, one has to diagonalize the energy operator in the space of 
admissible J. For F = 7/2 we have a 1 x 1 matrix, so the values given in 
Tables [KT] and [Bl2] for F = 7/2 are the energy levels shifts. The effect of the 
contact interactions is comparable with the uncertainty in the experimental 
value of the ^^N quadrupole moment, being two orders of the magnitude lower 
than the quadrupole-orbit interaction. 



C Natural widths and contact P-wave interactions in /i-meson ex- 
otic atoms 

Sections 3.1 - 3.4 describe interactions of nuclei and particles with arbitrary 
masses and spins. These results can be applied to /i-meson exotic atoms. In 
heavy nuclei, the Bohr radius in //-meson exotic atoms is smaller than the 
nuclear radii and the problem is relativistic in addition. We discuss therefore 
the finite volume effects and the relativistic effects affecting the natural widths 
and the contact P-wave interactions. 



C.l Natural widths of -meson exotic atoms in 2P states 

Due to the dipole 2P — 15* transition, the width of the 2P level equals 
rem 4aa;L , /'2\8 ^m'^ 
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where / = l-S", i — 2P, ujfi — Ef — Ei, and 




(C.2) 



It should be compared to the magnitude of the potential Uq^s ~ -^Q z'>Ti'^a{aZ)^ . 
Condition TJ^ < Uq^s gives roughly Qz > 10 /{aZm'^) - 6000/Z fm^. The 
highest electric quadrupole moments of nuclei are about 500 fm^, so it would 
make sense to check high-Z nuclei. 

The electromagnetic current has the form j = —eZpz/M + ep/m, where pz 
is momentum of the nucleus and p is momentum of the muon. In the center- 
of-mass frame, j = ep/m". The quantity ep/m' represents the convection 
current, which is the component of the total current ea, o: is the Dirac matrix. 
The nucleus spin current is neglected. We use expression j = em'cx./m" . For 
transition current, one has j/j = ieujfim'xfi/m" . The dipole transition matrix 
element x^j is calculated using wave functions obtained from solution of the 
Dirac equation. 

C.2 Electric charge and quadrupole moment densities in nuclei and nuclear 
electrostatic field 

Let us consider electrostatic potential, $o- created by the uniformly distributed 
electric charge inside of a sphere of radius R = 1.2 A^^^ fm: 




(C.3) 



so that we have for the electrostatic field 




r < R, 



(C.4) 



71°^, r>R. 



The charge density is given by 



47rpo = divEo = < 




r < R, 



(C.5) 



0, r>R, 



and normalized to 




(C.6) 
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The quadrupole component of the electrostatic potential has the form 



$2 = < 



|^ln(^)Qfn"n^ r < R, 



(C.7) 



i^gf n"n^, r>R, 



2r3 



where i?i = Re^^^. The radial dependence of $2 corresponds to the uniform 
radial distribution of the quadrupole component of the electric charge. 

The quadrupole electrostatic field E2 can be decomposed to the sum of L = 1 
and L — 3 components, the former constitutes the analogue of the delta- 
function component discussed in Sect. 3.4: 



-fln(f)Qfn^ r<R, 



(C.8) 



0, r>R, 
r < R, 



(C.9) 



5 (^PlrpaP-) ^ > p 

2^4 -'a ) ' ^ -n-, 



where 



The fact that the L — 1 component is localized inside the nucleus indicates 
that we deal with a contact interaction. The divergence of the L = 3 compo- 
nent of the quadrupole electrostatic field contributes to the charge density at 
r < R and the quadrupole moment also, so that we have 



47rp2i = divE2i = -^Qf n"n^, r < R, 



47rp23 — divE 



23 ■ 



15 



(C.IO) 
(C.ll) 



and P21 = P23 — for r > R. The normalization of $2 is chosen to satisfy 



(C.12) 



where p2 = P21 + P23- 
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C.3 Point size nucleus and nonrelativistic approximation 



The nonrelativistic reduction of the Dirac equation in external electrostatic 
field gives 

AU = -^^^^ ([E X p]cT - [p X E](T + divE) . (C.13) 

(for g = 2 see Ref. [65j). Substituting in this expression E2 = E21 + E23 and 
integrating the short-range part of the matrix element of AU at |x| < R one 
gets, in the limit of R ^ 0, 

= ^^^Qf (5r"/^(L, L) - r"^(L, S)X(O). (C.14) 

The numerical coefficients in front of the spin operators are distinct from those 
in Eq. (13.311) , since Eq. ( ]C.14I) includes the effect of the L = 3 component of E2 
and the Darwin term ~ divE. 



The radial part of the long-range component |x| > i? of the quadrupole inter- 
action diverges logarithmically at i? ^ 0. The spin matrix element, however, 
vanishes in L = 1 state, so AU^'^'^-range ^ q_ 

The Dirac equation contains the contact P-wave interactions. Relativistic wave 
equations for high-spin particles piB2] , obviously, contain the contact P-wave 
interactions also. The Darwin term for S = 3/2 has the form 

The coefficient can be restored using Eq.(9b) of Ref. [7]. 



C.4 Finite size nucleus and relativistic approximation 



The nonrelativistic limit of the Dirac equation ioi g = 2 and M — > 00 results 
to the quadrupole-spin contact P-wave interaction originating from the lower 
components of the Dirac bispinors. The lower components produce, however, 
other interactions too. We restrict ourselves with evaluation of the Larmor 
and Thomas contact P-wave interactions corresponding to E21. 

In order to discriminate the contact interactions, we use the Gordon's decom- 
position of the electromagnetic current: 

Jm = fr + = d, -2eA,) + |a,.(5 + A) ^, (C.16) 
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Table C.l 

The magnitudes of quadrupole-orbit long-range and quadrupole-spin contact P- 
wave interactions and natural decay widths of 2P states of muonic atoms formed 
with several spin I > 1 low- and high-Z nuclei. The experimental values of the 
nuclear electric quadrupole moments Qz are taken from Ref. [51], errors are not 

[2] 

displayed. UmLx is the maximum of the absolute value over F of the diagonal matrix 
elements < FJ\Uq^l\FJ > for J = 3/2 and n = 2, Umax is defined similarly for 
Ug^g. The effects of relativity and the finite volume of nuclei are included. Tgp is 
the radiation width of the 2P3/2 — 15'i/2 transition. 



Nuclei 


2H 


6Li 


^Li 


9Be 


10b 


iiB 


i*N 


/ 


1 


1 


3/2 


3/2 


3 


3/2 


1 


Qz [fm2] 


0.29 


-0.08 


-4.06 


5.3 


8.47 


4.07 


2.00 


ui^L [eV] 


4.4 X 10"* 


3.7 X 10"^ 


1.9 X 10-1 


5.8 X 10-1 


8.8 X 10-1 


8.8 X 10-1 


1.2 


UI^L [eV] 


5.8 X 10"^ 


4.2 X lO"'^ 


2.1 X 10-5 


1.2 X 10-* 


2.7 X 10-* 


2.8 X 10"* 


7.2 X 10-* 


[eV] 


8.8 X 10~^ 


7.3 X 10-3 


7.2 X 10-8 


2.3 X 10-2 


5.8 X 10-2 


5.8 X 10-2 


2.2 X 10-1 


Nuclei 


181 Ta 


185 Re 


190 Ir 


i98lr 


w^Au 


285 U 


253 Es 


/ 


7/2 


5/2 


4 


3/2 


3/2 


7/2 


7/2 


Qz [fm2] 


317 


218 


285 


75.1 


54.7 


493 


670 


uiL [oV] 


5.9 X 10-1 


5.5 X 10* 


5.3 X lO-i 


3.5 X 10-1 


2.7 X lO* 


1.2 X 105 


1.7 X 105 


Ull^ [eV] 


1.7 X 10^ 


1.7 X 108 


1.7 X 108 


1.1 X 108 


8.6 X 102 


4.6 X 108 


7.0 X 108 


r^™ [eV] 


7.2 X 10^ 


7.7 X 10^ 


8.4 X 10^ 


8.3 X 10^ 


9.0 X 102 


1.4 X 108 


1.9 X 108 



where = ($0, 0), $0 is given by Eq. (1C.3p . The time-like component of the 
spin current interacts with the quadrupole electrostatic potential: 



= / ^2jrdV. (C.17) 



The interaction energy flC.171) corresponding to constitutes the relativistic 
counterpart of the Larmor contact P-wave interaction discussed in Sect. 3.4. 
The relativistic extension of Eq. (13.291) . which takes the finite size of the nucleus 
into account, can be written as follows: 



U^^L,= .^J^^J .. r"^a, I)r"^(S, L) 



40mm' 2/(2/ - 1^ 



^ o , 25 , ,/?. 2m' 



X / r^rfr— In(-) .J hiir), 
W r e + m' — V 







where fujii^) is the upper radial component of the Dirac wave function in 
the potential V = e$o; with the normalization conventions of Ref. and e 
is the energy of the nJL level. In the limit of Z ^ and R —>■ 0, fujii^) ^ 
RnLir) and we recover Eq. (13.291) . In order to arrive at Eq. OC.lSp . we drop a 
term proportional to divE2i whose physical origin is attributed to the muon 
Zitterbewegung. It contributes to hyperfine structure of the P-wave levels and 
is included into Eq. (lC.14p . 



40 



The Larmor component of the interaction is locahzed completely inside of the 
nucleus. To get the correct isotope dependence in the limit of Z — >■ 0, one has 
to use the mass m in the electromagnetic current and the reduced mass m' in 
the Dirac equation. 

The interaction of the convection current with the quadrupole elec- 

trostatic potential generates the long-range quadrupole-orbit interaction and 
the Thomas component of the quadrupole- spin interaction. The latter comes 
from the lower components of the bispinors. The relativistic extension of the 
Thomas component of the quadrupole-spin contact interaction has the form 



R ^ 25 , .R^ Am'{e-V) 



/•« 2 r 25 , .K^ 4m' e- V ,2 , ^ 
Jo R^ £ + m'-\/ 2^"-^-^^ ^ 



Correct isotope dependence in the limit of Z — is reproduced by using the 
m' in the Dirac wave functions. The off-diagonal matrix elements receive an 
additional dependence on quantum numbers JL from the radial integrals. 

In fixed- J multiplets at Z — ^ and i? — 0, the ratio between strengths of 
the Larmor component of the contact P-wave quadrupole-spin energy and 
the long-range quadrupole-orbit energy (i.e. between the right-hand sides of 
Eqs.(B.20) and (B.18)) equals {aZ)'^/4. Due to the relativistic and finite vol- 
ume corrections, the width Ffp in ^ss-y j^gggives an additional factor of 0.22 
as compared to the nonrelativistic formula fIC.ip . the average value of 
entering the quadrupole-orbit interaction receives a factor of 0.42, the Larmor 
component of the contact P-wave interaction is suppressed by a factor of 0.15, 
and the Thomas component of the contact P-wave interaction is suppressed 
by a factor of 0.16. The numbers are given for 2P3/2 state. 

Results reported in Table [UH give an idea about the magnitudes of the natural 
widths and the contact P-wave interactions in muonic atoms formed with low- 
and high-Z nuclei. 



D Spin-1/2 and spin-3/2 relativistic spinors and SUq relations for 
octet and decuplet coupling constants 



The relativistic spinors of spin-1/2 and spin-3/2 particles are normalized by 



u{p,s)u{p,s) = l, 



(D.l) 
(D.2) 
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In the rest frame, these spinors have the form 



«=ri, (D.3) 







M° = I I , u 







(D.4) 



where (y? is the Pauh spinor and $, is the rest-frame spin-vector obeying the 
condition cr ■ ^ = needed to ehminate the spin-1/2 component from ^: 



^M = Eq:_.e(/i)e(--). (D.5) 



Here, the vector e*^™) with spin projection m is defined by 



e(-) = -^cT-^im), (D.6) 

where ^^^(+1) = ^^^(-1) = V2^^^i0) = ^2^2^(0) = 1, other components of 
spin-tensor ^°'^{m) vanish, (m) = Ca^y^^^im), and Ca/s = i^a"^)^. 

The spinors and spin-vectors with fixed spin projections are normahzed con- 
ventionally: 



^(^:(m')v5"(m) = 5^,^, (D.7) 

a 

ECK)rM = <5^,™. (D.8) 
The completeness conditions have the form 



2s + 



-Y.^-{m)^l{m) = -5l (D.9) 

^ m ^ 

^ I E r (m)a*(m) = 1 (515^^ - y^\<y^)i) . (D.IO) 



2s + 



Applying the boost transformation to the spinor indices of u and (see e.g. 
[65] . Chap. 3), and additionally, to the vector indices of using the matrix 



L(v) of Eq. (]A.4ll . one gets the relativistic spinors: 
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u{p, s) 
u^{p, s) 
u(p, s) 



lE + m I ^ 
2m 



(crp)y? 



lE + m I ^(pO 



2m 



lE + m i $, 

2m 



1 



E+m 



m{E+m) 



P(pl))y 



(D.li: 



(D.12) 



One can check that obeys p^u^{p, s) = and 7^u^(p, s) = 0. 



Expressions fID.lip and flD.12p can be used to get the nonrelativistic reduction 
Eqs.f HTTB]) - Km of the vertices NAK and ilEK. 

We use the SUq symmetry to fix sign of the interference term. Firstly, we get 
relations for the isovector pseudoscalar coupling constants gNNn and gANn and 
the isovector vector coupling constants Qnnp, Jnnp and f^Np- The relations 
for strange baryons can be obtained afterwards like in Sect. 4.3 with the use 
of the SU3 symmetry. 

The source of the pion field is divergence of the isovector axial vector current 

Jp(0) = apg"^<r3, (D.13) 

i 

where the summation extends to quarks, Ap is an unknown real constant, and 
is the pion momentum. 

In the A —>■ Np transition, the magnetic component is dominant. The source 
of the p mesons field of the magnetic type is rotor of the isovector axial vector 
current 

J^(0) = zAv^e-%^^a7r3, (D.14) 

i 

where the summation extends to quarks. Ay is an unknown real constant, and 
is the p-meson momentum. 

The matrix elements of the operator X]j (^T'^f over the proton and A"*" quark 
wave functions with the spin projections +1/2 are expressed in terms of the 
matrix elements of the corresponding baryons wave functions: 



1 15 



1 4 4? 

2 73^ V3 



(D.15) 
(D.16) 



43 



In the right-hand sides, a" act on the baryon spin indices. In deriving these 
equations, we used the A and proton wave functions constructed by compo- 
sition like in Eq. OD.Sp of the spin-1/2 d-quark and spin-1 uu-diquark wave 
functions. 

Comparison with the nonrelativistic matrix elements of Jp(0) and Jy(0) of 
Sect. 4 gives gNNn/C^rriN) = -|Ap, gANn/rriN = 2\/2Xp and {gNNp + Innp) 
I (2m at) = — fAy, m^f^^Np/^'N = 2v^Ay. In order to pass frompvr^ channel to 
p7r+ channel, we take into account the factor ^3/2; similarly for pp channels. 

Finally, we obtain 



JANp - 



3^2 



-gNNTT, 



rriN 



5 rriA 



{gNNp + Innp) 



(D.17) 
(D.18) 



The value of the AA^vr coupling constant is slightly away from the empirical 
value (14. 9 p for gNNn/V^^ = 3.67. The AA^p coupling constant is found to be 
fANp/V^ = -2.82 for gNNp/V^ = 0.84 and fNNp/V^ = 3.53 [M], model 
NCS97a, in the excellent agreement with Eq. (14.100 . 

The relative phases of the octet and decuplet coupling constants are thereby 
fixed. Coming back to Eq. (14.440 . we observe that the interference term TZi^ in 
the AS decay channel is positive. The remaining ambiguities affect phase of 
the total amplitudes, but not the interference. 

Assuming the SUq symmetry holds, the pseudoscalar meson exchnage and the 
magnetic vector meson exchange interfere in S-wave constructively: 

Consider first the nonstrange sector. In the nonrelativistic approximation, the 
PWBA amplitude of reaction AB —>■ CD due to tc^ and exchanges can be 
written as follows 



9Jt~< C|Jp(0)|A >< D\Jp{0)\B > 



< C\J^{0)\A >< D\J^{0)\B > 



m 



2 ■ 

V 



(D.19) 



Using Eqs. (1D.13P and (ID. 140 and averaging the amplitude over the directions 
of the momentum transferred q, one arrives at 



q2 + m*p 3 + rri' 
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As we mentioned, Ap and Ay are real constants (as a consequence of hermiticity 
of the currents), and so the two terms in brackets both are positive. 

If we replace rf by and by f/- and y-spin generators of the SU^ group, 
the statement on the constructive interference extends further to the charged 
vr- and p-mesons and kaons. 

The pseudoscalar and magnetic vector coupling constants are therefore propor- 
tional to the same quark operator. If the ratio between gNN-n and gNNp + Innp 
is taken positive, it remains positive for other members of the pseudoscalar 
and vector meson octets. The model NCS97a [53] fulfills such requirement for 
all coupling constants except for SSM and UNM, whereas the model NCS97f 
[53] fulfills it without exceptions. So, by following the model NCS97a we arrive 
at a destructive S-wave interference in the SS channel. The model NCS97f 
predicts a constructive interference there. The models NCS97a and NCS97f 
both predict constructive interference in the S'-wave dominant AS channel in 
agreement with SU^. 



E Angular part of gradient 



In the momentum representation, the angular part of gradient is defined 
as operator p-^ acting on functions of unit vectors n = p/p (see e.g. |33|, 
Chap. VII). The knowledge of identities listed below allows to simplify the 
calculation of integrals entering Eqs. fl4.55p - (14.581) : 



d 

dpa 

d d 
dpa dpp 



d d 



0, 



d 1 
op p 



1 d 
p op 

pZ 



2d I , , 

+ + ^ V V , 



dpa dpa dp"^ p dp p^ 



(E.l) 
(E.2) 
(E.3) 
(E.4) 

(E.5) 



(E.6) 

(E.7) 
(E.8) 



The integrals of angular variables, entering Eqs. (14.551) - (I4.58p . obey the fol- 
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lowing properties: 



J Yif^,Yimdnn = Sl>lSm'm, (E.9) 
/ (V">^W) iV'^Ylm) rffin = l{l + l)Sl'lSm'm, (E.IO) 

J {L-Y;^,) {L"YiJdQ^ = l{l + l)6in6^>m, (E.ll) 

/ (V" V'' Y;„,) (v" Ylm) = l\l + If5vi5m'mi (E.12) 

where I is the orbital quantum number. Equation flE.QP is the orthogonality 
condition for spherical harmonics Yjm(n), the next two equations give the 
normalization condition for the electric and magnetic spherical vectors. 
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